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In  this  papar  a  claaa  of  iterative 
procedures  la  presented  for  reconstructing  a 
finito  duration  sequence  from  no lay  aaaploa  of 
lta  Fourier  transform  phasa.  Thoaa  aaaauraatanta 
ara  ceablnad  with  a  knowladga  of  tba  truo  trans- 
fora  aagnituda  ahd/or  hard  constraints  on  tha 
phasa  noiaa  variations  to  dafina  sots  whosa 
intarsoction  must  contain  tha  trua  saguanco.  Tha 
algorithms  itarata  botwaan  tha  sequence  domain 
and  tha  ~  transform  domain  applying  tha  known 
constraints  (i.a.  finito  duration  and  known 
limits  on  phasa  variation)  in  oacb  domain. 
Results  of  an  experimental  investigation  are  pre¬ 
sented.  A  method  is  described  for  the  case  where 
limits  on  both  the  magnitude  and  phasa  variation 
of  a  finite  length  sequence  are  known. 


phase  information)  [4].  In  this  study  It  has 
bean  abown  that  if  x(n)  and  y(n)  ara  two  real, 
causal  (or  anticausal),  finite  duration  sequences 
and  if  certain  restrictions  on  tha  placement  of 
the  saros  of  X(s)  are  mat,  then  x(n)  is  equal  to 
y(n)  if  their  respective  signed  magnitudes  are 
equal  at  all  frequencies.  Unfortunately,  a  given 
finite  set  of  samples  of  the  signed  magnitude  is 
not  always  sufficient  to  uniquely  specify  a  se¬ 
quence.  This  is  true  even  through  it  is  known 
that  a  maximum  number  of  3H  suitably  chosen  fre¬ 
quency  samples  suffice  to  uniquely  specify  a 
sequence.  The  key  words  are  ‘suitably  chosen* 
since  the  required  smaples  vary  from  sequence  to 
sequence.  In  practice,  picking  a  sufficiently 
long  transform  length  (typically  10  times  the 
length  of  the  sequence)  allows  excellent 
restorations. 


In  recent  years  the  problem  of  signal  re¬ 
storation  from  partial  or  incomplete  information 
has  received  considerable  attention.  In  the  mag¬ 
nitude  retrieval  problem,  it  is  sssumad  that  the 
Fourier  transform  phase,  4(«)  *  argtx(u)]  ,  of  a 
finite  duration  sequence  is  known  exactly  at  B-1 
distinct  frequencies  in  the  interval  (0,v)  where 
M  is  the  known  extent  of  the  sequence .  under 
these  conditions  and  with  certain  restrictions  on 
the  placement  of  the  sefos  of  X(s) ,  it  has  been 
shown  [1]  that  iterative  procedures  exist  which 
will  converge  to  the  unique  solution.  These 
algorithms  Iterate  between  the  sequence  dcauiln 
and  the  transform  domain  applying  the  known 
constraints  (i.e.  finite  duration  and  known 
phase)  in  each  domain.  Convergence  has  been 
proven  -  within  the  framework  of  nonexpan live 
mapping  theory  (2). 

In  contrast,  recovering  a  finite  duration 
sequence  x(n)  from  a  knowledge  of  its  Fourier 
transform  magnitude  alone  has  proven  to  be  a  much 
more  formidable  problem  and  few  sequences  meet 
the  uniqueness  criterion  presently  available 
(3).  This  situation  has  prompted  a  number  of 
investigators  to  study  the  reconstruction  of  a 
sequence  from  samples  of  its  signed  magnitude 
(i.e.  the  transform  magnitude  and  one  bit  of 


When  the  phase  of  X(u)  la  not  known 
accurately  or  when  it  is  corrupted  by  noise,  it 
has  been  observed  that  the  magnitude  retrieval 
procedure  described  above  may  perform  badly. 
Other  techniques,  such  as  reconstruction 
averaging  [5]  or  minimum  cross-entropy  methods 
[6]  may  also  perform  poorly,  particularly  in  low 
signal  to  noise  ratio  environments.  It  seems 
likely  then,  that  additional  information  is 
required  if  we  are  to  achieve  more  acceptable 
performance. 

In  this  paper,  a  new  class  of  algorithms  is 
presented  which  achieve  significantly  better 
reconstructions  over  those  obtained  in  the 
past.  These  procedures  utilise  a  knowledge  of 
the  true  transform  magnitude  and/or  bounds  on  the 
phase  noise  variation  to  define  constraint  sets 
whose  intersection  oontalns  the  undistorted 
signal.  At  each  stage  of  the  Iteration,  the  cur¬ 
rent  estimate  of  x(n)  is  projected  (in  some 
fashion)  onto  the  oonstraint  sets.  In  all  cases 
to  be  examined,  the  finite  duration  requirement 
is  enforced  in  the  sequence  domain. 

Signal  Restoration  -  Roles  Free  Hwmo 

Let  x(n)  denote  a  finite  duration  one-dimen¬ 
sional  sequence  which  is  sero  outside  the  the 
interval  0<ndi-1.  The  s-transform  and  the 
Fourier  transform  of  x(n)  will  be  denoted  by  X(s) 
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and  X(w) .  Since  X(u)  la.  In  general,  a  complex- 
valued  function  of  w,  it  nay  be  written  in  terns 
of  its  real  and  imaginary  parts,  Xgfu)  and  Xj  ( u) 
respectively,  or  in  terns  of  its  magnitude  and 
phase,  IX («) I  and  r$  («)  respectively,  as  follows: 


X  (<■>)  -  XR(tt)  ♦  jXjIdl) 


!♦,<“> 


with:  ^(u)  -  tan-1  [Xj  :u)  AR(“)  ]. 

A  aet  of  conditions  under  which  x(n)  aay  be 
recovered  fraa  samples  of  its  Fourier  transform 
phase  is  contained  in  the  following  theorem  (1): 

Theorem:  let  x(n)  be  a  sequence  which 
is  aero  outside  the  interval  0<n«-1 
with  x(0)  *  0  and  which  has  a  (-trans¬ 
form  with  no  seros  on  the  unit  circle 
or  in  eonjugate  reciprocal  pairs. 

Then  samples  of  the  Fourier  transform 
phase  ♦jjlu)  at  N-1  distinct 
frequencies  in  the  interval  0<uK*  suf¬ 
fice  to  uniquely  (i.e.  to  within  a 
positive  constant)  define  x(n). 

A  special  case  of  this  theorem  arises  when 
the  phase  semples  of  X(u>)  are  equally  Spaced 
around  the  unit  circle  such  as  occurs  in  the 
Discrete  Fourier  Transform  (DFT).  Let  N  be  the 
length  of  the  DFT.  Then,  if  M>2N,  the  Iterative 
procedure,  which  replaces  the  estimates  in  each 
domain  by  their  known  values,  will  converge  to 
the  unique  solution  x(n)  for  any  initial  starting 
point  «0*r  [2].  It  la  this  algorithm  that  moti¬ 
vates  our  methods  for  dealing  with  noisy  phase 
measurements. 

dismal  Pastor  at  ion  -  ■oisr  Ihmmo 


When  noise  is  added  to  the  phase  measure¬ 
ments,  the  phase-only  iterative  procedure 
described  above  generally  produces  poor 
reconstructions.  Approaches  which  use 

restoration  averaging  IS]  or  minimum  cross 
entropy  methods  (6]  also  seem  to  give  poor 
results.  This  is  especially  true  of  these 
methods  in  low  signal  to  noiae  ratio  eases.  In 
this  section,  we  will  present  a  number  of 
iterative  schemes  useful  in  reconstructing  a 
sequence  given  noisy  phase  measurements.  These 
methods  have  been  found  to  give  significantly 
better  reconstructions  than  have  been  previously 
obtained  albeit  at  tbe  coat  of  higher  information 
requirements,  lore  the  added  cost  takes  the  fora 
of  knowledge  of  the  true  magnitude  and/or  bard 
constraints  an  the  phase  noise  variation.  As  in 
the  noise-free  ease,  the  iteration  proceeds  by 
projecting  the  current  estimate  of  x(n)  onto  tbe 
constraint  aets  in  both  the  sequence  domain  and 
the  transform  domain.  In  all  cases  to  be 
examined ,  the  constraint  aet  in  the  aequenoe 
domain  la  the  knowledge  of  the  first  (non-sero) 
point  and  the  finite  duration  requirement .  Tbe 
four  methods  to  be  described  differ  only  in  what 
oenetltstes  the  constraint  set  in  the  transform 


Given  that  u.  »  2vk/M,  k  «  0,1,  ...  ,M-1  are 
a  set  of  K  distinct  frequencies  with  M>2N,  each 
of  the  projection  operators  can  be  described  as 
follows  (see  Figure  1  for  illustrations  of  the 
constraint  sets  defined  below.): 

•  Method  I:  No  magnitude  information  is 
available  and  the  phase  noise  is  known 
to  vary  no  more  then  (  A  radians  from 
its  true  value.  Alternatively,  we  can 
say  that  the  true  phase  is  no  more 
than  s  A  radians  distant  from  the 
measured  ^oiay  value.  At  each  stage 
of  the  iteration,  we  examine  the  phase 
estimate  4(u^).  If  it  is  more  than 

A  radians  distant  from  the 
neSsurenent  value  we  replace  the  esti¬ 
mate  with  the  measurement.  The  magni¬ 
tude  is  not  altered. 

•  Method  II :  Bere,  the  true  transform 
magnitude  |X(u)|  is  known  for  all 
values  of  k,  but,  no  Information  is 
available  on  the  variability  of 

p(u^).  We  do,  however,  have  the 
noisy  phase  measurements.  For  each 
value  of  k  and  at  each  stage  of  the 
Iteration  we  replace  the  magnitude 
estimate  by  its  true  value.  Tbe  phase 
estimate,  however,  is  not  changed.  It 
should  be  pointed  out  that  this 
procedure  differs  from  the  phase 
retrieval  problem,  as  it  usually  takes 
form,  in  that  an  initial,  if  noisy, 
estimate  of  the  phase  is  available. 

•  Method  III:  Knowledge  of  the  true 

magnitude  is  combined  in  this  method 
with  a  hard  constraint  A  on  the  phase 
noise  variation  to  d&ine  an  arc 
centered  on  the  measured  phase 
value.  At  each  stage  of  the 
iteration,  the  current  estimate  of 
X<« )  is  projected,  using  a  nearest 
neighbor  rule,  onto  the  convex  hull  of 
these  arcs. 

O  Msthod  IV:  As  in  the  third  method, 
tbe  true  magnitude  |X(u  )|  is  known  as 
Is  the  maximum  phase  noise 
variation  A  .  However,  at  each  stage 
of  the  proves,  the  magnitude  estimate 
is  replaced  by  tbe  true  value  of  tbe 
magnitude.  Tbe  phase  estimate  is 
replaced  by  tbe  noisy  phase  measure¬ 
ment  if  the  estimate  lies  more 
than  A  radians  away  from  the  phase 
measurSnent. 

To  test  these  four  methods,  an  eight-point 
non  minimum  phase  sequence  (x(n)  -  (4,2, 

-11, 5, 4, 5, 15, -4) J  was  used.  An  M-point  DFT  of 
this  sequence  is  calculated  and  to  it's  phase  is 
added  an  M-point  noiae  sequence.  This  noise  se¬ 
quence  is  generated  using  a  uniform  probability 
density  function  given  by: 


-T7-  ,  -A  <W <A  , 

P6  («)  "  1  ^  P  P 
p  0  ,  otherwise. 

Each  of  the  four  Methods  were  exsained  for  trans¬ 
fora  lengths  of  M  «  16,  32,  64  and  128_  points  and 
noise  variations  A  equal  to  */2,  10"T»,  10"2*, 
10_3»,  10"**,  and  TO"5*.  A  relaxation  paraaeter 
(2]  of  0.99  is  used  in  these  experiaente.  Two 
perforaance  Measures  were  calculated.  They  are 
1)  the  no realised  aean  .  quare  error  (NMSE) 
defined  by  [5] 


l  [x(n)  -  Ax.(n)  ]2 

_ A  • 


l  x2(n) 
n-0 

and  2)  the  aean  square  error  iaproveaent  ratio 
(MSBI): 

M-1  , 

l  [x(n)  -  x  (n) J 
_  n-0 

MSE1  ■  - :  ■ 

I  [x(n)  -  Ax.(n)  ]2 
n-0 

In  these  expressions,  x(n)  is  the  true 
sequence,  x (n)  is  the  initial  or  starting  asti- 
aate  of  x(n),  xf(n)  is  the  final  estiaate  of 
x(n),  and  A  is  a  par  meter  chosen  to  ainiaise  the 
noraalised  aean  square  error. 

A  nuaber  of  observations  can  be  aade  froa 
the  results  obtained  thus  far  in  our  investi¬ 
gations: 

•  Methods  XI  and  XV  per  fora  phe.  beat  in  the 
noise  ranges  10"'*  to  10"*  s.  For  noise 
variations  greater  than  10"' i.  Method  IV  is 
the  algor ltha  of  choice.  For  exaaple,  if  A 
-  »/2,  Method  II  obtains  only  21  dB.  of  aean 
square  error  iaproveaent.  Method  IV,  on  the 
other  hand,  gives  a  nearly  perfect 
reconstruction  with-,  an  MSBI  of  242  dB.  for  N 
equal  to  64  points  and  5000  iterations 


•  At  intermediate  noise  levels  (i.e.  10"1*  to 
10"5*)  nethods  II  and  IV  give  roughly 
equivalent  results.  Bowever,  as  M  becoaes 
larger.  Method  TV  tends  to  eonverga  auch 
aore  quickly.  In  the  saae  range  of  noise 
levels.  Method  XI  and  XV  converge  at  the 
saae  rate  as  the  phase-only  iteration  using 
the  sequence  without  noise  added. 

0  Method  III  obtains  large  reductions  in  the 
error  in  the  early  stages  but  very  little 
iaproveaent  (if  any)  is  obtained  as  the 
Iteration  proceeds. 

0  Method  X  is  the  least  effective  of  the  four 
Methods  exaiined.  However,  it  also  requires 
the  least  aaount  of  a  priori  infornatlon 
about  x(n)  or  its  transfora. 

•  As  a  rule,  increasing  M,  the  length  of  the 
transfers,  laproves  the  oonvergenee  rate  of 


the  algor lthas.  Bowever,  in  soar  Instances, 
a  shorter  DFT  length  can  outperfora  a  longer 
length. 

Figure  2  Illustrates  these  observations 

concerning  the  four  algor lthas.  Here,  M  is  equal 
to  64,  and  for  those  sequences  with  noise  present 
it  is  given  a  aaxiaua  variation  of  10"  * 
radians.  Also  shown  in  this  figure  is  the 

perforaance  of  the  phase-only  Iteration  for  the 
phase  saaples  with  no  phase  noise  added  and  with 
a  noise  level  of  A  -  10"  ■  added.  The 

differences  between  tile  initial  aean  square 

errors  are  attributable  to  whether  or  not 
aagnitude  information  is  available. 

Shown  in  Figure  3,  is  the  perforaance  of 
Method  IV  as  A  is  varied.  81allar  effects  are 
seen  in  the  results  froa  the  other  procedures. 
In  this  figure,  M  is  equal  to  64. 


Beetoratloa 


In  this  section,  we  relax  the  need  for  know¬ 
ledge  of  the  true  aagnitude  and  require  instead 
that  bounds  on  the  variations  in  phase  and  aagni¬ 
tude  noise  be  known.  These  bounds  define  a 
region  about  the  Measured  aagnitude  and  phase 
which  aust  contain  the  true  transfora  values  of 
the  sequence.  At  each  stage  of  the  iteration  we 
project,  using  a  nearest  neighbor  rule,  the  esti¬ 
aate  of  X(w.  )  onto  the  convex  hull  of  the  region 
described  above.  Again,  in  the  sequence  domain 
we  enforce  the  finite  duration  requirement. 
Also,  we  assuae  x(0)  is  known. 


Figure  4  represents  the  results  of  an  ex- 
perlaent  using  the  eight-point  sequence  defined 
earlier,  A  equal  to  10"3»,  and  M  equal  to  64. 
Various  lfvels  of  uniformly  dlstr lbuted  noise 
ranging  from  10"'  down  to  10"5  were  added  to  the 
magnitude  data  and  reconstructions  were  obtained 
as  shown.  Also  shown  is  the  behavior  of  the 
phase-only  iteration  when  given  the  true 
sequence's  phase  and  when  given  the  true  phase 
corrupted  by  phase  noise  at  a  level  of  10"3*. 
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FIGURE  V.  CONSTRAINT  SETS  FOR  THE  FOUR  METH- 
OOS  DISCUSSED  IN  THE  PAPER:  la)  method  I  -  true  ae- 
quence  it  known  to  lie  in  the  intersection  of  M  wadga 
thaqad  regions,  (b)  method  II  -  true  taquanca  liat  in  the 
intersection  of  M  cirdae,  (c)  method  III  -  true  taquanca  Imm 
in  the  intartaction  of  the  convex  hulls  of  M  arcs,  and  (d) 
method  IV  —  true  taquanca  it  known  to  lia  in  the  intersect¬ 
ion  of  M  arcs. 
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FIGURE  2.  COMPARISON  OF  METHODS  (N-S.M-S4): 
(a)  phaaa  only  Iteration  -  no  noise,  (b)  phaaa  only  iter¬ 
ation  -  noiaa  level  of  0.001  tr ,  (c)  method  I  -  notes  level 
of  0.001  it ,  (d)  method  II  -  noise  level  of  0.001  ff  ,  (e) 
method  III  -  noise  level  of  0.001  if ,  (f)  method  IV  - 
noiaa  level  of  0.001* . 


FIGURE  3.  EFFECTS  OF  NOISE  LEVEL  ON  PERFOR¬ 
MANCE  USING  METHOD  IV  (N-8,M-64):  (a)  phase  only 
iteration  -  no  noise,  (b)  non*  level  of  0.1  ir  ,  (el  noise  level 
of  0.01  n  ,  (d)  noise  level  of  0.0017T  ,  (a)  noiaa  leva!  of 
0.0001s  ,  (f)  noise  level  of  0.00001  tr  . 


AND  MAGNITUDE  USING  CONVEX  HULL  PROJECTIONS: 
(a)  phaaa  only  procedure  -  no  noise,  (b)  phase  only  iter¬ 
ation  -  phase  noise  level  of  0.001  V,  (c)  magnitude  noiaa 
level  of  0.1,  (d)  magnitude  noise  level  of  0.001,  (e)  magni¬ 
tude  note*  level  of  0.00001,  and  (f)  no  magnitude  noise, 
(note:  unlew  otherwise  specified  there  is  a  phase  noise 


level  of  0.0011 r  .) 


